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Abstract 

We study the influence of the virtual pion cloud in nuclear matter at finite den- 
sities and temperatures on the structure of the p- and w-mesons. The in-matter 
spectral function of the pion is obtained within a selfconsistent scheme of coupled 
Dyson equations where the coupling to the nucleon and the A(1232)-isobar reso- 
nance is taken into account. The selfenergies are determined using a two-particle 
irreducible (2PI) truncation scheme (^-derivable approximation) supplemented by 
Migdal's short range correlations for the particle-hole excitations. The so obtained 
spectral function of the pion is then used to calculate the in-medium changes of 
the vector-meson spectral functions. With increasing density and temperature a 
strong interplay of both vector-meson modes is observed. The four-transversality of 
the polarisation tensors of the vector-mesons is achieved by a projector technique. 
The resulting spectral functions of both vector-mesons and, through vector domi- 
nance, the implications of our results on the dilepton spectra are studied in their 
dependence on density and temperature. 
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1 Introduction 



It is an interesting question how the behaviour of hadrons changes in a dense 
hadronic medium. One of the driving ideas is that chiral symmetry is ex- 
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pected to become partially restored with increasing density or temperature. 
Therefore the spectral functions of hadrons in their dependence on the ther- 
modynamic parameters of the hadronic environment are of special interest. 
They directly show the possible effects on the mass and decay properties 
of the particles. One of the experimental accesses to observe the spectral 
functions in-matter is provided through the study of electron-positron- or 
muon anti-muon pairs, called dileptons [1,2,3,4,5] both in hadron-nucleus 
and nucleus-nucleus collisions. The advantage of such electromagnetic probes 
instead of strong interacting particles, like pions or kaons, is that one gets 
information directly from the centre of the interaction zone. Above invari- 
ant masses of 400 MeV the pairs are mainly produced through the decay 
of vector-mesons like the p- and c<j-meson. Additional contributions come 
from bremsstrahlung and final-state Dalitz-decays of other resonances, like 
the 77-meson, which in part are controlled by other observations. Presently all 
nuclear collision experiments [2,3,4,5] observe a significantly enhanced dilep- 
ton yield in the invariant mass range below the light vector-meson masses 
compared to straight extrapolations from elementary processes, like proton- 
proton scattering. Several mechanisms, especially a lowering of the p-meson 
mass or an increase of its damping width, were proposed to explain the ex- 
perimental facts [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. Addi- 
tional experiments for photo-production off nuclei are planned, where one likes 
to investigate the u ir^'-f 37 channel [6,7]. Despite the perturbation 
through the final state interaction of the vr'^-meson [6] it will still be possible 
to isolate the a;-meson component because the competing p-meson process is 
highly suppressed. 

Ever since the early considerations of vector-mesons in dense matter [27,28,15] 
much effort has been devoted to study the properties of the p-meson, which 
strongly couples to the two-pion channel. However, with a few exceptions, e.g. 
[10], the in-medium properties of the pion are commonly neglected, just em- 
ploying unperturbed free pion states. Furthermore so far there are only few at- 
tempts to study the in-medium properties of the c<j-meson [17,18,19,20,21,22,29,30,31]. 
Thus, a more realistic determination of the in-matter pion cloud permits to 
investigate the dependence of the vector-meson properties on the various com- 
ponents of the pion spectral function in the medium. Compared to a schematic 
three level model for the pion modes studied by Wachs [31] we use a scheme 
of coupled Dyson equations in which the pionic degrees of freedom are mod- 
ified selfconsistently through the coupling to the nucleon and the A (1232)- 
resonance. All selfenergies are determined within a truncated two-particle ir- 
reducible (2PI) effective action formalism (^-functional) in their full depen- 
dence on energy and three momentum. This guarantees conservation laws of 
the theory to be fulfilled on the level of expectation values [32] despite the 
approximations employed. From the in-medium spectral function of the pion 
we are able to calculate the modified spectral functions for the p-meson as 
well as for the a;-meson in coupled scheme. Thereby the process p — > a; -|- tt 
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leads to a strong interplay between both vector-mesons. In particular the u- 
width is very sensitive to the space-like, i.e. low energy component of the pion 
caused by its coupling to nucleon nucleon-hole states. For a realistic treatment 
of this component Migdal's short range correlations [33] are important. These 
RPA-type correlations shift the spectral strength to higher energies, such that 
pion condensation is prevented. In order to concentrate on the effects due to 
the modified pion we neglected direct couplings of the vector-mesons to the 
baryonic currents. The latter have been investigated by several authors [8,9], 
who found a multicomponent structure for the in-medium p meson which is 
yet absent in our approach. 

Dealing with vector-mesons requires to consider current conservation on the 
correlator level. This is ensured by a projector formalism [34] which constructs 
a four-transversal polarisation tensor. 

The paper is organised in the following way: We first summarise the essential 
definitions in the realtime formalism used in the calculations at finite temper- 
ature and density. Then we present the model for the in-medium pion coupled 
to the the baryonic degrees of freedom. Afterwards we show how to derive a 
four-transversal and thus a current conserving polarisation tensor by means of 
a projector formalism. For the so obtained in-medium spectral functions of the 
vector-mesons we use the vector-dominance concept in order to calculate the 
effect of our results on the dilepton rates. The different decay contributions as 
well as the effect of possible mass-shifts of the vector-mesons in the dilepton 
spectra are investigated. A discussion of our results concludes the paper. 



2 Realtime formalism and 2PI generating functional 

Interested in all spectral properties as a function of real energies the realtime 
formalism [35] will be used rather than the Matsubara technique. Within this 
formalism all two-point functions such as the propagator ^ or selfenergies take 
four components, e.g. 

iG—{x, x') = (T^(x)^^(x')) iG++{x, x') = {T-^^{x)¥{x')) 
iG+-{x,x') = {^{x)'^\x')) iG-+{x,x') = t{^Hx')^{x)), 

depending on the order of the fermionic (upper sign) or bosonic (lower sign) 
field operators \t'(,T) appearing in the expectation values. The space-time coor- 
dinates are denoted by x and x', the time-, respectively, anti-time-ordering is 



^ We use units with c = h = 1 throughout the paper. 
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symolized by T and T ^. The Dyson equations then become matrix equations 
with respect to the {— , +} labels (placement on the closed time contour). 

For a uniform medium in thermal equilibrium one can exploit the standard 
simplification by going to the space-time Fourier-transformation x — x' ^ p 
of all quantities, where p = {p^,p) is the four momentum. Then the coupled 
Dyson equations simply become algebraic and it is sufficient to solve for the 
retarded (or advanced) Dyson equation, since all four components can be 
deduced from this by thermal Fermi-Dirac or Bose-Einstein factors arising 
from the Kubo-Martin-Schwinger (KMS) condition. Thereby the retarded and 
advanced Green functions follow from 

G^(p) ^ G-{p) - G-+(p) = G+-{p) - G++(p) 

G^ip) ^ G--{p) - G^-{p) = G-^{p) - G^+ip) (1) 

G^ip) ^ {G^{p)y 

with the KMS conditions 

iG-+ = TnT{p)A{p) 

iG+-^ {IT nT{p))A{p), ^ > 

where the spectral function and thermal weights are given by 

iA{p) = G^ip) - G«(p) = (G«(p))t - G«(p) 

1 (3) 



nrip) 



exp((pO-/x)/T)±l 



with temperature T, chemical potential // and the energy The above fisted 
relations for G likewise apply for all equilibrium two-point functions, i.e. also 
for the sclfcnergies. They are also written in a way that they apply to Dirac 
spinors or vector fields, whereby all quantities additionally become Dirac- 
matrix or Lorentz-tensor structures, respectively. The adjungation denoted 
by the f-symbol refers to this matrix structure ^ . For relativistic bosons the 
spectral functions A(p) are normalised to 27r with respect to the integration 
measure p° dp^. 

The propagators and self energies of the baryons are denoted by G and E, 
while we use D and U for the corresponding quantities for the meson fields. In 
selfconsistent two-particle irreducible (2PI) approximation schemes [32,36] the 
selfenergies S or 11 for the baryons and mesons are derived from a generating 
functional, called ^-functional. This functional is given by a truncated set 
of closed diagrams in accordance with the interaction Lagrangian where all 
fines denote dressed, i.e. selfconsistent propagators. The selfenergies result as 



^ Dirac spinors need to be treated with special care since one has to employ the 
"hermitian" form ^^ji^ of the Gamma-matrices before executing the adjungation. 
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functional variations with respect to the propagators, i.e. 



iT.{x,y) 



SiG{y,x) 



or 



iU{x,y) 



di^G,D] 
5iD{y,x) ' 



(4) 



This imphes an opening of a corresponding propagator hne in the diagrams 
of For the resulting set of coupled Dyson equations such 2PI approaches 
guarantee that even in a partial resummation of a single class of diagrams 
the conservation laws which are related to the symmetries of the system are 
fulfilled on the level of expectation values [32] . 



3 The Model 

The present investigations are restricted to spin as well as isospin symmetric 
nuclear matter. This implies vanishing chemical potentials for all mesons in 
the light flavour sector. All sclfconsistent calculations of the spectral functions 
are carried out numerically as a function of energy and momentum on a two 
dimensional grid of 5 MeV by 5 MeV/c resolution using an iterative procedure. 

3.1 Pions and baryonic Resonances 

In the baryon sector the nucleon and the A(1232)-resonance are included. The 
pion is considered as Goldstone boson of chiral symmetry which, in lowest ap- 
proximation in pion-energy, leads to p-wavc couplings to the baryonic pseudo- 
vector currents together with a minor s-wavc term which we ignore. For the 
energies involved we approximate the spinor structure of the baryons by the 
non-relativistic limit, while the kinematics is kept in the relativistic form. This 
approximation allows to handle the spin-3/2-structure of the A- resonance in 
a simple way ^ . Furthermore, for the investigated range in excitation energies 
we can safely use the no-sea approximation for the baryons. Thus we neglect 
processes involving anti-baryons and set the baryon spectral functions to zero 
for negative energies. 

The following non-relativistic form of the interaction Lagrangian 



(5) 



^ A fully relativistic formulation would require the use of special interactions as 
can be seen for example in ref. [37], which then would also have to be incorporated 
in the later introduced Migdal formalism. 
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is used with the spinor-field operators and of the nucleon and A- 
resonance. Here cr and r denote the vectors of Pauh matrices in spin and 
isospin, respectively. Likewise, the vector S defines the three spin-coupling 
matrices of spin 3/2 and 1/2 to spin 1, while T is the same in isospin-space. 
In our scheme we neglect pion self-interactions {^^^^ oc tt^) which give rise to 
selfenergy diagrams of sunset type [34] as well as contributions from AA~^ 
excitations. Both are of some importance at temperatures above 80 MeV, 
where they would even further broaden the pion spectral functions. 

Furthermore we supplement an exponential formfactor 



at the ttNN- and ttNA- vertices. Values of A= 440 MeV and ^f^ATA = 0.02 MeV"^ 
are required for a decent fit of the pion-nucleon phase-shifts in the 33-chanal ^ 
and a proper vacuum spectral function of the A-resonance. For Qt^nn we use 
0.007 MeV-^ [39]. 

For the chosen Lagrangian mean-field terms drop out and the first diagrams 
for the ^-functional in the baryon-pion sector which we use are: 



The lines represent the fully dressed propagators, while apart from the form- 
factor (6) the vertices remain bare. 

Not all of the selfenergy terms obtained through (4) are of equal importance. 
For instance the A-vr-loop contribution to the nucleon selfenergy is found 
negligible compared to the A^-7r-loop and thus neglected in the production 
runs. The retarded pion selfenergy can be obtained from a dispersion relation 
which in no-sea approximation does not require any renormalization 



Interested in the pion modes a complete determination of the baryon selfen- 
ergies E(p) is beyond the scope and not the aim of the present model. For 
simplicity we drop the real parts of T,{p), which would require a renormal- 
ization procedure within a selfconsistent resummation scheme [40]. However 
the normalisation of the corresponding spectral functions is restored in each 
iteration step. 

^ This implies a cut-off dependence of the in-medium results. In models with a 
larger body of resonances, of. [38], the cut-off A can be chosen above 700 MeV. 




(6) 
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As is well known, the coupling of the pion to nucleon-hole {NN~^) and AA^~^- 
states provides a strong softening of the pion modes already at normal nuclear 
densities. Since there are no indications from experiment for a pion condensate 
at normal nuclear density one has to include a repulsive force which shifts the 
spectral strength up to higher energies. This can be achieved by a repulsive 
short-range interaction between the baryons first introduced by Migdal [33] 

XXX 

For the Migdal parameters we used gu = 0.5, 5122 = 0.6 gi2 = 0.5 without 
any form factor. At a momentum scale of 200 MeV/c for the relevant NN~'^ 
excitations''' this choice is compatible with the recent parameter set of Suzuki 
[42] and Nakano [43] who calibrated it on Gamow- Teller transition strengths. 
These correlations are to be included in the retarded pion selfenergy in form 
of an RPA-resummation with particle-hole loops 



""^--00--^ -OOO-— 

where also mixtures with the A-resonance, i.e AA'^~''^-loops need to be con- 
sidered. In the non-relativistic limit for the spinor structure of the baryons, 
only the spatial components of the Dirac matrices survive in all loops which, 
by symmetry arguments, are diagonal (for a fully relativistic treatment cf. 
[45]). Therefore, apart from a q^-factor due to the difference between the 
pion-baryon coupling and the Migdal interaction the basic loops for the RPA- 
resummation are given by the normal pion selfenergy diagrams resulting from 
the ^-functional (4) 



ttR 

^ V NN- 



-O- -O- 



Thus, the RPA-resummation can be done without further problems. With 
the three possible vertices one finally obtains [38,11] 

= J 7Z-2 T-PFi? T7—2 —Rr ^ -jiTm TTr — ^ l^J \P) 



Ah 



for the retarded pion selfenergy. The corresponding results together with those 
for the pion spectral function are shown in the Figs. 1 to 4. 



^ These excitations, which essentially influence the p-LO mixing in matter, are mainly 
determined by gu. 
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Fig. 1. Pion spectral function A^^ at T=0 
MeV and p = po for two different mo- 
menta 
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Fig. 2. Real (dotted line) and imaginary 
(full line) part of the pion selfenergy at 
T=0 MeV and p = po for |p| = 200 MeV 
(upper panel) and |p| = 400 MeV (lower 
panel) 
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Fig. 3. Same as Fig. 1 at T=120 MeV Fig. 4. Same as Fig. 2 at T=120 MeV 



In the selfenergy and spectral function we see the effect of the coupling to 
the different l:)aryonic excitations. The on-shell pole broadens and is shifted 
towards lower energies due to the net-attraction of the interaction in this 
kinematical region. This component strongly mixes with the AN^^ component 
which leads to the shoulders around this peak. This component appears in the 
time-like region for low momenta, while it traverses to the space-like region 
for higher momenta. At low energies an entirely space-like pion component 
appears due to NN~^ excitations. From the kinematics it is clear that this 
component stays in the space-like region for all momenta. Such space-like pion 
modes have to be interpreted as a scattering process on the baryons in the 
matter mediated by the exchange of pions. Due to this component, the pion 
spectral function has non-vanishing strength at all energies, such that in all 
processes where pions are involved all thresholds disappear. Due to the short- 
range Migdal repulsion there is no spectral component where the peak position 
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decreases towards lower energies when one increases the momentum. The form 
of the spectral function turns out to be very sensitive on values for the Migdal 
parameters and the formfactor F. In addition, in the literature one finds quite 
different procedures to incorporate the formfactor in the Migdal resummation. 
Some of them include formfactors also for the Migdal- vertices [41,42,43] ^ . We 
followed the approach from rcf. [11] with formfactor-free Migdal vertices. The 
difference between these two techniques turns out to have a great effect on 
the shape of the pion spectral function for momenta in the order of 400 MeV. 
For higher temperatures all structures are smoothed out and the value of the 
selfenergy even decreases because of limited number of available nucleons to 
scatter. 



3.2 Pions and Vector-mesons 



In this chapter we study the influence of the medium-modified pion-modes on 
the spectral functions of the p- and c<j-mesons. The widths of these resonances 
are generated through the decay into two pions in the case of the p-meson, 
respectively three pions for the a;-meson. For the a;-meson we choose the in- 
direct decay via the p-meson (the so called Gell-Mann, Sharp and Wagner 
(GSW) process [47]) because, as shown by Theileis [46], this is the dominant 
contribution in vacuum. The interaction Lagrangian defining the vertices is 
given by [48,49,50] 

^iS; = gpnnP^iT^d^^) + ^.p.e"^'^"a;„o'^p^a,7r, (9) 

for a review see [51]. In this sector we are lead to the following ^-functional 
approximation 




From the point of view of a ^-derivable approximation we neglected the cou- 
pling of the p-a;-loop back to the pion. This process is expected to contribute 

only marginally to the low energy region of the pion spectrum, which turns out 
to be the most relevant part to determine the width of the a;-meson. Therefore 
only the following retarded selfenergies are included: 



^ Following the Korpa-Malfliet [41,44] prescription for the Migdal vertices we could 
qualitatively reproduce their results as a check for our numerics. 
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R 



I 



\ 



n 



(11) 



/ 



\ 



n 



R 

7/1/ OJ/OTT 




(12) 



The correlation between p- and a;-meson modes results from the a;-meson 
selfenergy (12) and the reverse process encoded in the coupling of the p-meson 
to the cu-TT-loop 



which both are included selfconsistently. Since the model omits higher lying 
degrees of freedom the in- medium changes of the real parts of the vector-meson 
selfenergies are less precisely determined within the model space. Therefore 
we drop them, while restoring the normahsation of the spectral functions at 
each iteration step in the further analysis. In a later section wc discuss the 
sensitivity of the results on a possible change of the vector-meson masses in 
matter. 

3.3 Transversality of the vector-meson polarisation tensors 

The resulting polarisation-tensors (11 - 13) have to be four-transversal because 
of current conservation. In normal perturbation theory this is guaranteed or- 
der by order. In the Dyson approach, where one sums up a restricted subclass 
of diagrams to infinite order, one generally violates Ward identities on the 
correlator level. Thus the polarisation tensor may contain four-longitudinal 
components n['^(g), which have to be absent and which may lead to the prop- 
agation of unphysical degrees of freedom. From general grounds, this deficiency 
can be cured by corresponding vertex corrections. Without further approxima- 
tions, however, this leads to a presently intractable scheme of Bethe-Salpeter 
equations which accounts for the required i-channel exchanges required by 
crossing symmetry. We circumvent this problem in the following way. From 
transport considerations [52] it is known that such polarisation tensors have at 
least two relaxation times. Because of charge conservation, one of these times 
has to be infinite, implying that the component n°'^(g) vanishes exactly for 




(13) 
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q = and qo ^ 0, while the second relaxation time is clearly finite. Such a re- 
sult can never be reached in a truncated Dyson-resummation scheme where all 
relaxation times are finite, because they are determined by the damping-time 
scale of the dressed propagators involved in the loops. On the other hand, the 
spatial components of the polarisation tensor, given by the autocorrelation 
of spatial currents, have solely finite and short correlation times which can 
be expected to be safely approximated within a Dyson resummation scheme. 
Therefore our strategy assumes the spatial components of the polarisation ten- 
sors n^j^ to be given by the selfconsistent loops, while the time-components 
are to be corrected such that the full tensor becomes four-transversal. This is 
achieved by using a projection technique where the full tensors (11 - 13) are 
decomposed into a four-longitudinal part 11; and two four-transversal parts 
and Ily which are three-longitudinal and three-transversal, respectively. 



n'^'^(g) = nr (g) + nf (?) + (g) 



q^ J ^ ' (14) 



Tl^%q)=U^-^YTlT{q). 



Here g^^ is the metric tensor, whereas 5^^ and q^ have vanishing time compo- 
nents. We do not specify the IT; part any further, since this part has just to be 
dropped due to current conservation. The scalar functions 11^ and 11^ can be 
calculated solely from the spatial parts of the polarisation tensors using the 
following traces 



3n, = TV 



q^ q^ 

^-gikIl"'^2UT+^^-Ur. or 



n 



q^ 



(15) 



Ht = - (SHs - Hi) . 



Due to the traces, the complicated Lorentz-structure of the interaction ver- 
tices (9) in the loops can be reduced entirely to expressions involving the 
four-, respectively, three- vector squares q^,p^, p'^ and q^, etc. of the external 
and internal momenta involved in the loops. Furthermore, three-longitudinal 
and transverse modes decouple such that the corresponding spectral functions 
Ai^q) and v4r(g) are directly given by the scalar functions and 11^(5) 

of the corresponding components of the polarisation tensors. 
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3-4 Vector-meson spectral functions in matter 



With the dressed pion spectral function from section 3.1, the three-longitudinal 
respectively three-transversal polarisation tensors and spectral function of 
both vector-mesons can now be calculated. In Figs. 5 to 12 we show the results 
for normal nuclear density at different temperatures^ 




200 400 600 800 1000 
inv. mass [MeV] 

Fig. 5. 3-longitudinaI /9-meson spectral 
function Ap^i at p=200 MeV and p = po 
for different temperatures 
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/9-meson selfenergy \mYi^^ at p=200 
MeV and p = po for different temper- 
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Fig. 7. 3-transversal p-meson spectral 
function ^p,T at p=200 MeV and p = po 
for different temperatures 
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Fig. 8. 3-transversal part of the p-meson 
selfenergy Imll^-jn at p=200 MeV and 
p = Po for different temperatures 



^ Please note that by definition, cf. (14), Imll^ vanishes on the light cone 
and becomes positive in the space-like region. 
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Fig. 9. A^^L at p=200 and p = po 
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Fig. 10. Imn^^^ at p=200 MeV and 

P = Po 
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Fig. 11. A^^T at p=200 MeV and p 
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Fig. 12. Imn^,^ at p=200 MeV and 
P = Po 
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The main effect towards finite densities and/or finite temperatures is the dis- 
appearance of thresholds present in the free particle kinematics for the pi- 
ons in the loops. The spectral strength starts right at zero energy for both 
vector-mesons, though less visible for the a;-meson due to the smallness of 
the coupling. This has significant consequences for the low mass region of the 
corresponding dilepton yields (cf. next section). The temperature has no great 
infiuence on the p-meson spectrum, while the tu-meson width increases signifi- 
cantly. For the latter effect indirectly the nucleon-nucleon-hole excitations are 
of special importance as discussed later. 
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Fig. 13. e+e--rate at T=60 MeV p = po 
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Fig. 15. /x+/x~-rate at T=60 MeV p = po 
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Fig. 14. e+e"-rate at T=120 MeV 
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Fig. 16. p+p'-isXe at T=120 MeV 



3.5 Dileptons 



Using vector-dominance we can calculate [53,54] the production-rate for dilep- 
ton pairs (e"'"e~ and fi'^ 



dR a 



2 



d'^qd'^xdt 6n^q'^ \ 



q2 _ 2777,2 



(16) 

V i/vTrTrV ^''^yvTrTT/ 



4 



from the in-medium spectral functions of the p- and cu-meson. Here nB{qo,T) 
is the thermal Bose- Einstein weight, while niy and gfyTrir denote the vector- 
meson mass and the coupling constant of the meson to two pions. The mass of 
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the lepton, electron or muon, is symbolised by m/. In our model we obtain the 
results displayed in Figs. 13 to 16 for the direct decay of a vector-resonance 
into a dilepton pair via an intermediate time-like photon. 

With applications to nuclear collisions in mind, we use the vacuum expressions 
for the photon- and lepton-propagators, since these particles interact only 
weakly, such that they do not get modified by the surrounding medium. In 
addition, we used the coupling constants from [51] in order to describe the 
partial decay width of both mesons into dileptons. One can sec that, due to 
the thermal weight and the broad spectral functions, the dilepton spectra are 
essentially faUing over the whole energy range. This would clearly be different 
if the vacuum spectral function for the o^-meson had been used here. The 
/x"'"/i~-spectrum has a threshold at twice the muon mass. At higher invariant 
mass both di-lepton spectra merge, because then the muon mass becomes 
negligible in comparison with the energy. The contributions of the a;-meson 
are hardly visible in the total spectra because of its increased width. For equal 
widths of both mesons the contributions of the a;-meson falls below that of 
the p-meson by one order of magnitude because of the smaller partial decay 
width into dileptons. 

In the following we explicitly analyse the influence of the various components 
of the pion spectral function on the vector-mesons and thus on the dilepton 
spectra. Formally we do this by splitting the spectral function of the vector- 
meson into the various components related to the different processes feeding 
into this vector-meson channel. Thus, decomposing the total damping width 
into partial widths Fv^totlp) — Si ^v,i{p) the dilepton yield can be brought into 
a Breit-Wigner like form with partial in- and out-widths 



dR 3 

nTiqo)J2M<l) (-2Imn,, 



d'^qd^xdt (27r) 



4 



3 / N ■^-^ 4 q'qFv i Fy 



, (17) 



(27r)^ ^^^^ fr {q' - mlY + qlTl,,, 



(suppressing the tensor structure of spectral function and polarisation tensor 
which leads to a degeneracy factor 3 for vector particles). Here, Fy^e+e- is the 
dilepton decay width of vector-meson v. 

We start the discussion of the effects on the a;-meson. In the medium there are 
three major processes contributing to its damping width, illustrated by per- 

turbative "time-flow" diagrams where the time is running from left to right 
and vertical lines denote a virtual space-like propagator, which mediates a 
two-body interaction: 
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Fig. 17. e~^e~-rate from the decay of the 
w-meson at T=120 MeV p = po divided 
into different contributions: 

F: 7rp — > a; 
D: p ^ nu 
S: pN ujN 
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Fig. 18. e~''e~-rate from the decay of the 
p-meson at T=120 MeV p = po divided 
into different contributions: 

A: TTTT — > p F: ttcu — > p 
B: TT ^ 7rp D: a; — > 7rp 
S^: ttN ^ pN S^: uN ^ pN 



(F) = ssmj (D) = 

Here the subsequent decay of the CL-'-mcson into the virtual time-like photon 
and its final decay into the lepton pair is not illustrated. In the selfconsistent 
calculation all these processes are included automatically by using dressed 
propagators. For the first process (F), pvr — a; — > e"'"e~, the a;-meson is formed 
by the fusion of a p-meson with a quasi-real, in-medium pion. Its inverse exists 
already in vacuum and determines the vacuum decay width of the cj-meson. 
The second process (D), p — * ttcj — ^ 7re"'"e~, corresponds to a p-Dalitz-decay 
via an intermediate c<j-meson. In the selfconsistent calculations, both above 
mentioned processes just differ in the sign of the pion energy in the Trp-loop 
of the a;-selfenergy (12). The process (S) in (18) corresponds to the scattering 
pN— > a;N mediated by a virtual, i.e. space-like pion exchange. In view of 
the pion modes at zero temperature (cf. Fig. 1) we isolate this space-like 
component by a cut on the far space-like region with pion loop momenta with 
\p\ > 2|p°| in (12). At T = 120 MeV this separation is less evident in view of 
the broad structure of the pion spectral function (Fig. 3). Thus the different 
components of the processes displayed in Fig. 17 somewhat depend on this 
cut. In addition to the fusion width (F), which constitutes just a temperature 
dependent modification of the vacuum width (e.g. accounted for by Schneider 
and Weise [29]) a genuine in-medium process, namely the scattering process 
(S), contributes with comparable strength at the nominal resonance position. 



(S) 



(18) 
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The "Daliz" -decay of the p- meson (D) dominates the low mass region. 

In summary of this analysis: a major portion to the uj spectrum results from 
processes (S) which are not accounted for in the simple on-shell treatment. 
These contributions, however, sensitively depend on the in- medium properties 
of the virtual pion cloud, which certainly needs further clarifying investigations 
before quantitative conclusions can be drawn. 

The same type of processes as in (18) also occur for the p-meson just inter- 
changing p with cj, listed in the p-meson decomposition given in Fig. 18 as 
(F), (D) and (S^^). However, these c<j-induced components are less important 
compared to the coupling to the two-pion channels 




At invariant masses above 300 MeV, the 7r+7r~-annihilation process (A) is 
clearly dominant. Processes (B) and (5*-^) are not present in any on-shell 
treatment of the pion as they solely arise from genuine off-shell components 
of the pion spectral function. In this respect (B) can be interpreted as a 
bremsstrahlung process radiated off a pion scattered in the medium, while 
process (S) corresponds to inelastic ttN pN scatterings mediated via vir- 
tual pion exchange. The latter two components, which emerge completely 
consistently within the model, only contribute to the very low-mass part of 
the invariant mass spectrum. 

3.6 Dependence on density and temperature 

We restrict the discussion of the dependence of the vector-meson properties on 
density and temperature to the case where the baryons and the vector-mesons 
retain their vacuum masses and use a cut-off A = 440 MeV (6) at the pion- 
nucleon- and pion- nucleon- A- vertices (5). The results for the damping width 
at resonance condition, i.e. at the vacuum masses, are summarised in Figs. 19 
to 22. 

Both vector-mesons show a strong non-linear dependence on the baryon den- 
sity p induced by the selfconsistent treatment. For the p- meson the relative 

changes are less dramatic due to the large vacuum width of 150 MeV from 
the TTTT-annihilation channel. The temperature effects are different. Here, with 
increasing temperature, the p-meson width stays nearly constant at nuclear 
saturation density po. The main contribution to the increase of the total width 
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Fig. 20. at T=30 MeV depending on 
the density. 
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comes from the processes p un and pn — > u. On the other hand the uj- 
meson is seen to have a strong dependence on temperature inducing a width 
change of up to a factor of 3 in the temperature range from T=0 MeV to 
T=120 MeV. For both mesons the temperature effects are due to an enhanced 
scattering with the pions in the medium . In comparison to the pcrtubative 
calculations by Schneider and Weise [29] the selfconsistent treatment produces 
an even broader cu spectral function, due to the appearance of the low energy, 
space-like, pion modes. The QCD sum rules approach of Klingl and Weise 
[56] and also the low-density expansion scheme of Lutz et.al. [14] produced 
cj-mcson widths of about 40 MeV at T = and p = po, which is nearly 2/3 
of ours. The slight deviations to our results can be understood by the change 
of the vector-meson masses in medium in their models as discussed in chapter 
3.7. Compared to the (j-a;-mixing model of Saito et al. [30] we cannot account 



^0 This has already been pointed out by Roy et al. [21] who also indicated that the 
process a; — > p7r might be increased in-matter due to mass shifts. 
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for level-level repulsion in our model, as we do not include the real parts of 
the vector-meson selfenergies. 

Altogether the a;-meson width sensitively depends on the properties on the in- 
medium pion cloud. An increase of the cut-off scale A of the TrA/^A/^-formfactor 
(6) for instance would significantly enhance the RPA pion-modes and induce 
a further increase of the cj-meson width. Once the cu-meson width increases 
beyond 50 MeV, it loses its prominent peak-structure and its component can 
no longer be easily resolved in the dilepton spectrum due to the dominance of 
the p contribution. 



3.7 Dependence on the vector-meson masses 



Our description of the vector-meson sector is not as complete as to permit 
definite predictions for the mass shifts of the vector-mesons, i.e. the real parts 
of the vector-meson selfenergies. This is a subject of current debate and in- 
vestigations and requires a more complete model space within hadronic de- 
scriptions [8,14] or support from QCD, e.g. in form of sum-rule constraints 
[22,23,24,25,26]. For simplicity the real parts of the vector-meson selfenergies 
are put to zero in the results presented above. On the other hand we can ex- 
plore the effect of mass shifts by simply changing the masses by hand. For the 
following discussion we concentrate on the damping widths (longitudinal and 
transverse) of the vector-mesons at the in-medium mass position = for 
a typical momentum of 200 MeV/c given by 

rv,L/T=-imnjL/T(?)/go . (20) 



A phenomenological lowering of the p-meson mass, e.g. within the Brown-Rho 
scaling scheme [12], has been considered by many authors. This leads to an 
enhanced strength in the mass region below the p-meson vacuum mass. From 
the calculations thus far the mass dependence of the vector-meson damping 
widths at T = and saturation density can approximately be read off by 
i.e. inspecting Figs. 8 and 12. They predict marginal changes for the p-meson 
width. Thus, the entire p-meson strength will be shifted, e.g. towards lower 
masses, with the extra benefit that the thermal weight factor in (17) enhances 
the dilepton yield by about a factor exp(Amp/T) . The a;-meson shows sensitive 
changes of its damping width as a function of its invariant mass. Since the only 
channel, to which the a;-meson couples to in our model, is the 7r-p-channel, we 
expect to depend sensitively only on the mass difference or ratio between 
the p- and a;-meson masses. 

The explicit results for the dependence on m^/mp for T = and normal 
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Fig. 26. e+e -spectra resulting 
from the decay of the cj-meson at 
= 0.85ma,,^ac., T=120 MeV and 
p = po, For notations see Fig. 17. 



nuclear density can be found in the Figs. 23 and 24 for the a;-meson and for 
the partial width Tp^^^^ of the p-meson, respectively. Due to the dominance 
of the TTTT-annihilation channel, the latter has minor impact on the total p- 
meson width. The displayed dependence follows expectations from phase-space 
arguments. Our results compare favourably with the calculations of Klingl and 
Weise [56] and Lutz et.al. [14], who used QCD sum rules or the low density 
theorem for the determination of the mass shifts. Although in the Klingl and 
Weise result [56] the mass of the a;-meson is lowered by 15% relative to the 
p-meson mass, which would infer a much smaller width, one has to take into 
account that in those calculations a much broader p-meson (Fp ^ 300 MeV) 
was used inducing a compensating increase of the cu- width. We also studied the 
vector-meson mass dependence at finite temperature. In Fig. 25 we compare 
the width of the a;-meson with a mass reduced by 15% to the width without 
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such a reduction . We find that the overall offset is nearly the same for 
all temperatures. According to this, the relative effect decreases for higher 
temperatures. This can be understood because Idnematical effects, like mass 
shifts, are of less importance for broad spectral functions. In Fig. 26, we show 
the effect of this mass shift on the dilepton spectra. We observe that the peak 
structure becomes more prominent because of the smaller width of the oo- 
meson and that the total rate is enhanced. This enhancement is due to the 
thermal weight which increases exponentially towards lower temperature. 



4 Conclusions 

We have investigated the in-medium effects on the light vector-mesons due 
to the modification of the pion modes in nuclear matter at finite density and 
temperature. To isolate this effect only the coupling of the vector-mesons to 
the pion modes has been considered so far. The direct coupling to baryonic 
currents along the lines of refs. [8,9,14] within a selfconsistent scheme will be 
the subject of a forthcoming study. 

Within a selfconsistent Dyson resummation scheme, the nucleon and A(1232)- 
resonance were included as the main degrees of freedom in the baryon sector. 
Besides the direct ttNN and ttA^A couplings, accounted for in selfconsistent 
selfenergies up to one-loop order, also the short-range RPA correlations of 
Migdal type [33] were included in order to get a realistic behaviour of the low 
energy pion modes. While for the meson-baryon couplings a non-relativistic 
approximation was employed, the kinematics of all particles was treated rela- 
tivistically. The so obtained in-medium spectral functions of the pion provide 
the source for the vector-meson selfenergies through the Tnrp and Trpuj cou- 
pling vertices. While the former already induces the strong damping width 
of the p-meson in vacuum, the latter becomes essential for the strong vector- 
meson "mixing" at finite densities and temperatures. Consequences for the 
resulting in-matter dilepton spectra from the electromagnetic decay of the 
vector-mesons have been discussed. 

The polarisation-tensor for the vector-mesons has to be four-transversal such 
that no unphysical degrees of freedom are propagated. This was achieved by 

a projector technique (see Eq. (14) ff.) which provides the three-longitudinal 
and three-transversal parts, respectively. These two parts are identical in the 
vacuum due to Lorentz invariance, while they may differ in the medium. 

As an important result of our investigations we find a strong mixing of both 
vector-meson modes in-matter. It is induced by the low-lying space-like pion 

We show the mean value of longitudinal and transverse widths Tl and Tt- 
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modes which essentially arise from the coupling to particle-holc excitations. As 
a consequence, the cj-meson damping width drastically broadens both, with 
increasing baryon density and temperature reaching values of 100 MeV or even 
above. 

The quantitative results on the in-medium properties of the vector mesons 
still depend on the model assumptions and parameters for the in-medium 
pion physics such as the Migdal parameters and the formfactors for the nNN 
and ttA^A vertices. Despite more than two decades of work in this field, the in- 
medium properties of the pion are not yet conclusively settled. Thus, further 
improvements, especially concerning the A'"A'""^-component, are required for a 
quantitative understanding of the in-medium broadening of the light vector- 
meson spectral functions, especially of the c<j-meson. Here other techniques like 
the recent calculations from Lutz and Korpa [55] are to be mentioned which 
may offer alternative and quantitative strategies to understand the spectral 
function of the pion and thus provide a reliable basis for the calculations of the 
vector-mesons properties. As our approach qualitatively incorporates most of 
the relevant contributions to the pion spectral function, all qualitative features 
should already be visible in our results. Even if the situation at zero tempera- 
ture and saturation density may eventually become quantitatively settled, the 
extrapolation of the effective interactions to finite temperatures may requires 
further adaptions. 

The broad cj-meson distributions, as found in our calculations, will make it 
difficult to isolate this in-medium component in dilepton spectra. It will be 
detected in competition with the much stronger electromagnetic decay rates 
of the p-meson. Therefore complementary experiments which directly isolate 
the a;-component e.g. partly through hadronic channels, such as in the recent 
TAPS-experiments [6,7], are vital to clarify the situation. In nuclear collision 
experiments, one further has to face the fact that contributions from final- 
state interaction always arise. Such asymptotic state vector-mesons show the 
signature of the corresponding vacuum spectral functions and thus may hide 
the here addressed in-matter components. 
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